function h=errorbarfga(x, y,e)

if nargin < 3
  e = y;
  y = x;
  x = 1:length(y);
end

h=patch(...
   [x(1:end-1);...
    x(1:end-1);...
    x(2:end);...
    x(2:end)],...
   [y(1:end-1)-e(1:end-1); ...
    y(1:end-1)+e(1:end-1);...
    y(2:end  )+e(2:end);...
    y(2:end  )-e(2:end)],...
          [0.5,0.5,0.5],'EdgeColor','none'); 

hold on

plot(x, y, 'k', 'LineWidth', 2)
